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Abstract 

We study D-branes transverse to an abelian orbifold C^/.Z'„ x Z^. The moduh 
space of the gauge theory on the D-branes is analyzed by combinatorial calculation 
based on toric geometry. It is shown that the calculation is related to a problem to 
count the number of ground states of an antiferromagnetic Ising model. The lattice 
on which the Ising model is defined is a triangular one defined on the McKay quiver 
of the orbifold. 
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1 Introduction 



The moduli space of the gauge theory on D-branes transverse to an orbifold C^/F has been 
investigated based on the procedure proposed by Douglas, Greene and Morrison [1, 2, 3, 
4, 5]. In this approach the geometry of the moduli space, which is defined by F-flatness 
and D-flatness conditions, is analyzed in the language of toric variety and the result is 
described by a cone in a integer lattice Z'^ . The dimension c is determined by a calculation 
of a dual cone which is necessary to make the F-flatness and D-flatness conditions on the 
same footing. To find the dual cone we must carry out a combinatorial algorithm based 
on the linear programming, but it is difficult to obtain an explicit result in general. By 
inspecting closely the calculation of the dual cone, however, one can see that it does not 
fully use the structure of the gauge theory. 

One of the notable feature of the gauge theory is that it is a quiver gauge theory [6, 7]. 
Information on the gauge theory such as F-fiatness and D-fiatness conditions are encoded 
in the McKay quiver of the group F [8, 9]. So if we fully make use of the information as 
a quiver gauge theory, the calculation of the dual cone will become easier. Based on this 
observation we show that the combinatorial calculation of the dual cone is equivalent to a 
problem of counting the number of ground states of a certain Ising model on a lattice defined 
on the McKay quiver of the orbifold. In this paper we consider the orbifold C^/ Z„ x Z„, in 
which case the Ising model is an antiferromagnetic one on a triangular lattice with certain 
boundary conditions. 

In section 2 we review the gauge theory on one D3-brane transverse to a quotient 
singularity C^/F with F an abehan subgroup of SU{?>). We also recapitulate the procedure 
[1] for the analysis of the moduh space M. of the gauge theory. In section 3, we investigate 
the orbifold jZn x Z„. We show that the calculation of the dual cone can be replaced 
by a calculation of ground states of an antiferromagnetic Ising model on a triangular lattice 
constructed from the McKay quiver. In section 4, we present a systematic way to count 
the number of the ground states by using transfer matrices of the Ising model. We also 
generahze the analysis to the orbifold C^/Z„ x Z^. 

This paper is based on a talk "A relation between quiver variety and Ising model" at 
the autumn meeting of the Physical Society of Japan, September 2003. 

2 The moduli space for one D-brane transverse to an 
abelian quotient singularity C^/F 

In this section, we briefly review the gauge theory on one D3-brane transverse to a quotient 
singularity C^/F with F an abelian subgroup of SU{?>). 

2.1 The gauge theory on the D-brane and the McKay quiver 

The gauge theory on one D3-brane transverse to a quotient singularity C'^/F is obtained by 
considering |F| D3-branes on C^, and then projecting to C'^/F. The worldvolume theory on 
|F| D3-branes on is a theory with C/(|F|) gauge symmetry, containing a vector multiplet 
Aa (a = 0, • • • ,3) and three chiral multiplets X*^ (// = 1, 2, 3) in the adjoint representation 
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of the gauge group. The components of the gauge field Aa that survive the projection are 
those which satisfy 

^{g)Aa^ig)-' = Aa (2.1) 

where g & T, and ''y{g) is the regular representation of F. For an abelian group F, the 
regular representations takes the form 

lig) = ©rio'^.(^?) (2.2) 

where are one-dimensional irreducible representations of F. The condition (2.1) then 
implies that all elements of Aa must be zero except for diagonal components and hence 
the gauge group is isomorphic with [/(1)I'"L Note that the overall U{1) acts trivially on all 
fields, so the effective gauge symmetry is C/(l)l^l~^. The surviving components of are 
those which satisfy 

^{g)X'^^ig)-' = RigXX-' (2.3) 

where R{g) is the three-dimensional representation defining the action of F on C^. We 
choose the three-dimensional representation of the form 

R{g) = R,, (g) © R,, (g) © R,, (g) (2.4) 

with i?^^ © R^2 ® to be the trivial representation Rq. All elements of must be zero 
except for the 3|F| components X^ X^ X^ which we denote as 

Charges of the v-th U{1) gauge symmetry of Xu),yw, are given by 

?t)('^ti)) — ^W,V-\-V\ ^W,Vi 

respectively. 

Gauge group and matter content of the world volume theory on the D-brane are specified 
by a McKay quiver. The McKay quiver is a graph with a set of nodes V and a set of 
arrows A connecting nodes. Each node, which has one to one correspondence with the one- 
dimensional irreducible representation i?^ of F, gives rise to a factor t/(l) of the gauge group. 
The set of arrows of the quiver is determined by the following irreducible decomposition 

R0 R^ = ®nvu,Rw (2.7) 

Non-negative integer n^y^ is the number of arrows from node v to node w. For the repre- 
sentation R{g) in (2.4), the decomposition takes the form 

R(S)R,, = R,,+^, (g) © R^+,,, (g) © R^+y, (g), (2.8) 

thus we have three arrows (a„, 6„, c^) leaving the node Rv'. 

bv '■ Rv ^ ■Rv+V2 1 (2-9) 

. R^ ^ Ry^y^ . 

Note that arrows ay,by, c„ correspond to the fields Xy,yy,Zy, respectively. 
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2.2 The moduli space of the gauge theory 

The classical moduli space M. of the gauge theory is obtained by imposing the F-flatness 
and D-fiatness conditions of the theory and further dividing by the gauge group C/(l)l'"l~^. 
The F-flatness conditions come from the superpotential 

W ^Tr[X^,X^]X^ (2.10) 

whose minimization yields the matrix form of the F-term equations 

[X'',X'']^0. (2.11) 

In components, the equation takes the form, 

•^vVv+vi yv'^v+V2i yv^v+V2 ^vVv+v^i ^v-^v+v^ •^v^v+v\' (^■■^^) 

Note that the number of independent F-term equations is 2|r|— 2, and hence these equations 
define a |r| +2 dimensional complex subvariety Z. Since Z has the structure of an algebraic 
variety defined by a collection of monomial relations, it is an affine toric variety and admits 
an alternative presentation as a holomorphic quotient. If we change notation temporarily 
as 

(ui, • • • , liain) = (xo, • • • , X|r|-i, yo, • • • , y|r|-i, ^o, • • • , ^|r|-i), (2.13) 

we can solve the equations and express ui in terms of |r| -\- 2 parameters Va { a — 
I,--- ,|r|+2); 

a 

The rows of the 3|r| x (|r| -|- 2) matrix K = (Kia) are vectors in the lattice M = Z'^'"*"^, 
and define the edges of a cone a. To describe Z as a holomorphic quotient, we define a 
dual cone a as follows; 

= {n G RI^I+^ I m ■ n > for Vm e a}. (2.15) 

We introduce a matrix T such that its columns span the dual cone a. The dimension of the 
matrix T is (|r| -|- 2) x c where c is the number of edges of a. Each column of T corresponds 
to a homogeneous coordinate px {X = 1, ■ ■ ■ ,c) oi the holomorphic quotient description of 

Z = C7(C*)^-I^I-'. (2.16) 

The action of (C*)'^~'l'"l~^ on pa is determined by a (c — |r| — 2) x c matrix Qp defined by 
Qf = (KerT)*. 

Note that the variables are solved in terms of p^; 

Va^llpl"'. (2.17) 

A 

Combining these equations with (2.14), we obtain the relationship between ui and px 

n = llPx"'^"'- (2.18) 

X 
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Due to the definition of tlic dual cone, tlie exponent KiaTax = {KT)i\ takes value in Z>q. 

To obtain the vacuum moduli space M. of the quiver gauge theory, one must further 
impose | T | — 1 D-fiatness conditions that exist in the quiver gauge theory from the beginning. 
The D-term equation for the v-th U{1) group is given by 

3|r| 

?^(«OI«/r = = 0, • • • , |r| - 1) (2.19) 

1=1 

where C,^ is the Fayet-Iliopoulos parameter satisfying C,^ — 0. In the graph-theoretic 
language, the |r| x 3|r| matrix d with entries d^^i — qv{ui) is the incidence matrix of the 
quiver, where d^^i is —1 if v is the tail of /, +1 if v is the head of / and zero otherwise. 

Since the gauge group is f/(l)l^l~^ as noted above, one of the D-flatness conditions in 
(2.19) is redundant. So if we define a (|r| — 1) x 3|r| matrix A by deleting the last row of 
d, the |r| — 1 independent D-fiatness conditions are written as 

3|r| 

Y,^v,i\ui\^ = Cv. (^; = 1,--- ,|r|-l) (2.20) 
1=1 

To combine these conditions to the expression of Z (2.16), we must rewrite these condi- 
tions to those with respect to the coordinates pa, that is, we must find the U (l)'^'"^ charges 
for p\. The charge matrix Qu for px is determined from the charge matrix A for ui and 
the relation between the coordinates ui and p\ (2.18) as 

Qd{KTY = A. (2.21) 

The total number of D-fiatness conditions to be imposed on the auxiliary gauge theory 
is (c-|r|-2) + (|r|-l) = c-3, and gauge symmetry is U {ly-^^^-^ x U {ly^^-^ = U{iy-^. 
Thus the moduli space M takes the form 

M = {(pi, ...,Pc) e C^lic - 3) D-flatness conditions}/[/(l)'=-^ (2.22) 

Note that the |r| — 1 D-flatness conditions in equation (2.19) have Fayet-Ihopoulos param- 
eters 9i, while c — |r| — 2 D-flatness conditions coming from F-flatness conditions do not 
have such parameters. 

Finally the moduli space can be obtained in the form: 

M^(Z- S^)/(C*f^-^ = (C^ - ZJ/(C*)^-^ (2.23) 

where 5*^, are certain exceptional sets. The action of (C*)'^~^ on C is summarized in 
the (c — 3) X c matrix 

- ( S ) • f^-^^* 

and the columns of the 3 x c matrix Gt = Kcr((5t)* play the role of toric generators of A4. 
Before ending this section, we show that the matrix Gt is written as 

Gt = UKT (2.25) 
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if we introduce the 3 x c matrix 

/I ■■■ 1 ■■■ •■■ 
n= 0--- Ol--- lO--- 0|, (2.26) 
\0 ••• ••• 1 ••• 

where each row has |r| one and 2|r| zero. Since the matrix (Gt)* is defined as a kernel of 
Qti what we have to show is 

Firstly, the equation Qpi^lKTy = QpT^K^Tl* = follows from the equation QpT* = 0, 
which is obtained by the definition oi Qf- Secondly, 

QD{T^KTy = QD{KTfli* - AH* 

\ Q\T\-i{xw) E«, Q\r\-i{yw) E«; Q\T\-l{Zw) ) 
since the U{1) charge of x^, for example, satisfy 



0. (2.28) 



|r|-i |r|-i 

w=0 w=0 

Thus the matrix Gt takes the form 

/ m(x^) ■ n(pi) • • • m(,T„) ■ n{p^) \ 
Gt = IlKT=\ E.m{y,)-n{p,) ■■■ Ev m{y,) ■ n{pc) (2.30) 
\ E^, fn(z^) ■ ■^(Pi) • • • E^ 'fn(zy) ■ n(pc) J 

and the generators of the toric diagram for Ai are written as 

/ E^m{x^) -nipx) \ 
fix=[ E.m(|/„)-n(pA) . (A = l,---,c) (2.31) 

What wc would like to stress is that to obtain n\ we only need the inner product between 
m{zy) and n{;p\) and the expression n\ itself is not necessary. Based on this observation, 
we propose a method to calculate the matrix Of without solving the dual cone. Our strategy 
is as follows; we start with the following relation between u and p 



u^^\{pl^\ (2.32) 



A=l 



although the number c is not determined yet. Next we impose the F-flatness condition 
on the indices aix e Z>q and count all possible combinations of indices. The number c 
is determined as the number of independent combinations which satisfies the F-fiatness 
conditions. 
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3 CV^n X 

In this section we consider the gauge theory on the worldvolume of a D3-brane transverse 
to the quotient singularity C^/Z„ x Z„ [10, 11, 12, 13]. We show that the calculation of 
the dual cone is replaced by a tiling problem of a torus by rhombus, and a certain dimer 
problem. (Similar obserbation was made in [14].) We furthermore show that it is equivalent 
to a calculation of ground states of an antiferromagnetic Ising model on a triangular lattice 
defined on the McKay quiver. 



3.1 Configuration of indices on the McKay quiver 

The group x has one dimensional irreducible representations and apart from 
the overall C^(l), we have a [/(I)" ~^ gauge symmetry. If we choose the three-dimensional 
representation asi? = i?io®-Roi® -Rn-in-i) the generators g'1,5'2 of the group Z„ x Z„ 
act on as 



g2 : {X\X\X^)^{X\ujX\uj-^X^), 

where uj — . Then we have the following decomposition 

R ® Rij = Ri+i j ffi Ri j+i ffi Ri-i j-i, 
and hence there are three types of arrows 



(3.33) 



(3.34) 



Qij 


Rij 


Ri+lj-i 


bij 


Rij 


Rij+li 


Cij 


Rij 





(3.35) 



The McKay quiver of the quotient singularity C^/Z^ x Z^ is depicted in Figure 1. 



^00 ^12 ^10 

o — 



«oiO- 



o- 



^0 



•o- 



•o- 



^22 ^00 

■O 



►O ►O 



o 




O ►O 



o 



Figure 1: McKay quiver of the quotient singularity / Z2 x Z2 



This model has 3n^ variables Xij,yij,Zij, which corresponds to the arrows Qij, bij, Ci 
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respectively. The F-flatness conditions on these variables take the form 

•^ij Vi+lj Vij -^ij+ly 

Vij ^ij+l — ^ijVi-lj-li (3.36) 

As explained in the last section, we introduce new variables px {X = 1, ■ ■ ■ ,c) which are 
related to the original variables Xij, yij, Zij = 0, • • • , n — 1) by the following equations, 





a}- (x^ ■ 
= Pi" V2' ■ 


a?. 

yc , 




Vij 


= Vi V2 ■ 


■■ Pc\ 


(3.37) 


Zij 




7=, 
Pc , 





where the indices o;^-, P^j, 7^^ e Z>o are subject to the F-fiatness conditions. If we consider 
the quantity Xij i/i+ij -Zj+i j+i, we obtain the following relations 

•^ij Vi+l j Zi-\-\ j-\-\ Uij Xi jj^\ Zi^\ jj^\ 

— Vij Zijj^iXi^ij (3.38) 

Zij Vi—lj — l j ) 

due to the F-fiatness conditions (3.36). In terms of the indices cXij, PijjJij, these relations 
are written as 

'^ij + Pi+l j + 7i+l j+l — Pij + '^ij+l + 7iH-l j+1 

= + + (3.39) 

— lij + Pi-lj-l + j- 

To investigate sets of indices a^j, f3^j, 7^^- satisfying these conditions, it is useful to consider 
configurations of indices on the McKay quiver. As noted above, each variables Xij, Ui+ij 
and Zi^ij+i corresponds to arrows a^-, and q+i^+i. Thus in terms of the arrows, the 
relations (3.38) are written as 

O'ij bi+ij Cj_|_i = bij cti j+i Cj+i j+i 

= bij Ci j+i aj_i j (3.40) 

— Cij bi—ij—i (li—ij. 

These relations imply that all closed paths encircling triangles of the McKay quiver are 
equivalent as depicted in Figure 2. 

If we assign the indices afj, (5^j and jlj to arrows aij, j and Q+ij+i respectively, we 
obtain a configuration of indices on the McKay quiver. The conditions (3.39) imply that 
the sum of indices on every triangle of the McKay quiver must be equal. Thus we only 
need to consider configurations of indices satisfying this condition. We denote the sum of 
the indices of each triangle by s (See Figure 3). 

8 



o ►o ►o 




o ►o 



Figure 2: F-flatness conditions imply the equivalence among triangles in the McKay quiver. 



1 1 1 1 2 1 

1 1 1 1 

(^1^3^ ^1^3^ 
(a) (b) 

Figure 3: Configurations of indices on the McKay quiver for x Z'2. (a) is a config- 

uration with s = 3. (6) is not allowed since sums of indices are different from triangle to 
triangle. 



r^yr'T^ r^^'T^ VTr^y^ v^'^ 

11001 00000 01000 10001 

01100 00000 01000 00100 

^1^3^ ^1^1^ C^O^l^ C^O^l^ 

(a) (b) (c) (d) 

Figure 4: A configuration with s = 3 is written as a superposition of three configurations 
with s = 1. 
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Here an important point is that any configuration with the sum s can be written as 
a "superposition" of s configurations with s = 1 (See Figure 4). Based on this fact we 
assert that each configuration with s = 1 corresponds to a coordinate p\ and hence c is the 
number ofdifferent configurations with s = 1. Thus what we have to do is to enumerate all 
configurations of indices with s = 1. 

As an example we depict all configurations of indices with s = 1 for C^/Z2 x Z2 in 
Figure 5. 




Pi Pi P3 Pa. Ps Pi Pi P% P9 



Figure 5: Configurations with s = 1 on the McKay quiver for C^/Z2 x Z2. Bold lines 
represent links assigned index one and thin lines represent links assigned index zero. Each 
configuration corresponds to a coordinate px. 



There are nine configurations with s ~ 1, which means that c — 9. By using the 
equation (3.37), we obtain the relation between Xij, i/ij, Zij and pi, • • • ,P9 as follows; 



xoo = PiPePs, xio = P1P7P9, xoi = P1P6P9, 2^11 = PiPvPs, 

y00=P2P4P8, yiO=P2P4P9, Voi ^ P2P5PQ: Vu ^ P2P5P8: 
Zoo = P3P4P6, Zio = P3P4P7: ^01 = PsPsPg, ^11 = PsPbPl- 



From the equation (2.18), we obtain a matrix KT as 



KT 



(mi ■ nx 



f 




n2 




n4 


ng 


Ue 


n-j 


ns 


ng 


mi [xoo) 


1 














1 





1 





m2 (xio) 


1 

















1 





1 


{xoi) 


1 














1 








1 


rui {xii) 


1 

















1 


1 





"^5 (yoo) 





1 





1 











1 





(yio) 





1 





1 














1 


^7 (l/Ol) 





1 








1 











1 


"^8 (yii) 





1 








1 








1 





"^9 (^00 ) 








1 


1 





1 



















1 


1 








1 








mil {zoi) 








1 





1 


1 











\ mi2 {zii) 








1 





1 





1 









This result is the same as that obtained by the calculation of dual cone. 
By defining the charge projection matrix 11 as 

111100000000 

n= 1000011100000 

000000001111 



(3.41) 



(3.42) 



(3.43) 
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we obtain nine toric generators rix according to the equation (2.30), 



G, = UKT 





hi 


h2 


"-3 


n4 


ns 


ne 


hj 


ng 


ng 


\ 




4 














2 


2 


2 


2 









4 





2 


2 








2 


2 




V 








4 


2 


2 


2 


2 








/ 



(3.44) 



Thus we finally obtain the toric diagram as depicted in Figure 6. 



P2 

P4 Ps * P9 



P3 P6 Pi Pi 



Figure 6: Toric diagram for the model x Z2. We assign a coordinate p\ to each 

toric generator. 



3.2 Enumeration of configuration of indices on the McKay quiver 

Now we would like to present a systematic way to count the configurations of indices with 
s = 1. As we will see below, the configuration of indices is translated to other systems such 
as a tiling problem, a dimer problem and an Ising model. 

Firstly, we show that enumerating configurations of indices with s = 1 on the McKay 
quiver for C^/ x Z„ is equivalent to a certain tiling problem. What is tiled is a torus 
with size n x n and tiles are three types of rhombuses depicted in Figure 7. 




X - type y - type z - type 



Figure 7: Three types of rhombuses. Each rhombus should be put on the McKay quiver 
such that its diagonal lies on the indices with one. 



Each rhombus should be put on the McKay quiver such that its diagonal lies on the link 
with indices one and its boundary on the link with indices zero. Then the conditions (3.39) 
on the indices assures that rhombuses cover the torus and do not overlap each other. The 
fact that configurations of indices with the sum s is a superposition of s configurations with 
s = 1 corresponds to the fact that a s-tiple tiling is decomposed to s tilings with s = 1. As 
an example, we depict all tilings for the model C^/ Z2 x Z2 in Figure 8. 
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Pi Pi Pi Pa Ps Pi Pi Pi P9 

Figure 8: Tiling a torus with size 2 x 2 by rhombus. 



Secondly, the configuration of indices is translated to a dimer problem on a hexagonal 
lattice. The vertices of the lattice lie on the center of the triangles of the McKay quiver 
and dimers connect neighboring vertices such that they intersect to links with indices one. 

Thirdly, the problem to enumerate configurations of indices with s = 1 is equivalent 
to a problem to count ground states of an Ising model on a lattice defined on the quiver 
diagram. Lattice sites of this model are the nodes Rij of the McKay quiver, and spin 
variables ± at the nodes interact through neighboring links connected by arrows. Thus the 
lattice is a triangular one with n x n sites on a torus. We consider an antiferromagnetic 
Ising model on this lattice. Although an antiferromagnetic Ising model favors anti-parallel 

links (H — ) and ( — h), there must exist parallel links (++) and ( ) due to frustration 

since the lattice is triangular. Therefore every triangle in the lattice has two parallel links 
and one anti-parallel link in the ground states. By deleting parallel links from the lattice, 
we obtain a diagram of anti-parallel links. Anti-parallel links form rhombuses because of 
the property of the ground state of the Ising model. Thus each ground state of the Ising 
model corresponds to a tiling by rhombuses discussed above. 

Note that not only periodic boundary condition but also ant-periodic conditions are 
allowed. Thus there are four boundary conditions 

CTj+nj; (^ij — (^ij+ni (3.45) 

where "p" and "a" represent periodic and anti-periodic boundary conditions. An impor- 
tant point is that interchange of -|- spin and — spin gives the identical configuration, the 
correspondence between ground states of the Ising model and tilings is two to one. 

Ground states of the Ising model on the McKay quiver for C^/ Z2 x Z2 is depicted in 
Figure 9. The boundary conditions for {^1,^2,^3}, {^4,^5}, {Vq^Vi} and {ps^Pg} are (p,p), 
(p, a), (a, p) and (a, a), respectively. 

4 Counting the number of the ground states of the 
Ising model 

In this section we count the number of ground states of the Ising model introduced in the 
last section. 



m ZZ 



(p, p) Gij = 

(p, a) aij = 

(a, p) aij = 

(a, a) (Tij = 
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Pi Pi Pi Pa 



Ps Pi Pi Pi P9 



Figure 9: Ground states with (Too = + of the antiferromagnetic Ising model on the lattice 
of the McKay quiver for C^/ Z2 x Z2. 



The energy E of the Ising model on the triangular lattice is given by 

n-l 

E ^ J ^ (o-yO-j+y + aijaij+i + ai+ij+iaij). (4.1) 

i,j=0 

For an antiferromagnetic Ising model, the coupling constant J is positive. The partition 
function of the Ising model is 

Zn.n = E (4.2) 

where the sum is taken over all configurations of spins which satisfy (anti-) peri- 

odic boundary conditions in the horizontal and vertical directions. Therefore the partition 
function is decomposed into a sum of partition functions with the four boundary conditions, 

7 _ 7(PP) I 7(pa) _|_ 7(ap) I yiaa.) (a q\ 

In Z^^l, for example, summation is restricted to configurations of spins {aij} with the (p, 
p) boundary condition. 

To calculate Z^^l, we introduce a transfer matrix T^^^ as follows; 

n— 1 n—l 

{<7ij} j=0 i=0 

n-l 

= J2U(Ti'%<^... (4-4) 

Wij} j=o 

= Tr(T^P))". 

The suffix aj of the transfer matrix Tn^^ represents a set of spin variables of the j-th row 
{(Toj, • • • , an-ij}, so T^^^ has 2" x 2" elements. We use the lexicographic order for the suffix 
of the matrix. Using the transfer matrix T^^^ , the partition function for the (pa) boundary 
condition is written as 

4':<t = Tr[(T(P))"i?„], (4.5) 

where the 2" x 2"^ matrix 
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is inserted in order to realize the anti-periodic boundary conditions in the vertical direction: 

(Tjn ^ -CTjO- (4.7) 

Similarly, the partition functions for the boundary conditions (ap) and (aa) are written as 

zifn = IV(^(-))^ (4.8) 

= Tr[{T^tYRn], (4.9) 

(a) 

where T„ is the transfer matrix for the anti-periodic boundary condition in the horizontal 
direction: 

O-nj = -O-Qj- (4.10) 

What we have to do next is to find the transfer matrices T^f^ and Ti^\ For this purpose, 
we consider configurations of spins on the 2 x n lattice depicted in Figure 10 (a) with free 
boundary conditions. 



(a) 



;+i ^2 .,+1 



6 - 

<^0 J J J 



j + l ^n-l j+l 



Q 



— O 



(b) 



^0 , + l ,+1 ^2 j+l 



— 6 - 

^'o J J ^2 j 



o- 



— o 

^n-2J ^n-lj ^Oj 



(c) 



CJq ; + 1 f^l , + 1 ^2 j+l 



^0 j ^1 j ^2 j 



^n-2 j+l ^n-1 j+l ^0 j+l 

- — 



o 



-o 



^n-2j ^n-lj -%j 



Figure 10: lattice 



The 2"^ X 2"^ transfer matrix r„ corresponding to this configuration is written as 

n-2 
i=0 



This matrix satisfies a following recursion relation 

{TnU, = {Tn-l H^),,,{Hn-2 ® T2)c 



(4.11) 
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where the 2*^ x 2^^ matrix if^ is given by 



Hk — 



1 1 
1 1 



By using the exphcit form of the matrix T2: 



To 



with K — pj, we are able to obtain the matrix T„. 



th^^ tfc^^ 







1 


1 




\ 




















e^' 


e-2^ 




v 


e^ 


1 


1 


e-3i^ 


/ 



(+hH+) 




(4.12) 



(4.13) 






C+>H--) 




tB^^ tB^^ 




t^^i) (z-y-^j t^^j t^^. 



C+>^) C+>h9 (^^^^) 




1 1 e"'*^ 

Figure 11: The matrix elements of 72- 



Configurations of spins with the periodic boundary condition in the horizontal direction 
are realized by adding a pair of spins uoo and aoi on the right side of the lattice as depicted in 
Figure 10 (b). Thus the transfer matrix Tn^^ is obtained from T„ by adding the contribution 
from the last plaquet in Figure 10 (b). Therefore we obtain 

with 

{Pn)cr a- = e~^'^[^^'^"~^-'''^°^"'"'^"~2^'+l'^°^'+l^"'"'^"~2^''^"~^^'+l"'"^°^'+l^"~^^] (4 15) 

By using the matrix T2, the 2" x 2" matrix P„ is written as 

P„ = Xl{T2 (8) Hn-2)Xn (4.16) 
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with 



x„ = (X2®ir^)(i2®x„_i), X2 



For example, the matrix T2'^^ is given by 

/ 



/I 
















1 








1 








Vo 








1/ 



(4.17) 



.(p) 





1 


1 




1 




^2K 


1 


1 




^2K 


1 


^2k 


1 


1 





/ 



(4.18) 



The transfer matrix tIi''' for the antiperiodic boundary condition afong the horizontal 
direction is obtained from by adding the contribution from the last plaquet in Figure 10 
(c). Thus Tn^ is written as 



with 



(A \ _ ^—l3j[-^{—Crn-2jO-Oj—(7n-2j+lO-Oj+l)+Cr„-2jan-2j+l—Croj+lTn-2j^ 



The 2" X 2" matrix A„ is given by 



where 



l2"-l 

I2T1-1 



(4.19) 
(4.20) 

(4.21) 
(4.22) 



makes the boundary condition antiperiodic in the vertical direction. For example, the 



(a) 

matrix T2 is given by 



Tr 



(a) 



^-2k 


^2K 


^-2K 


^2K 




^-2K 


^2K 


^2K 


^2K 


^2K 


^^2K 


^-2K 


g2fc 


g-2K 


^2K 


g-2K 



(4.23) 



Combining these results, the partition function is written as 

Z„xn = 'IV[{(TiP))" + (Ti^))n(l + Rn)\. 

For example, we obtain 



^2x2 = 18e^^ + 32 + 12e-^^ + 2e-^2^. 



(4.24) 



(4.25) 



The first term proportional to corresponds to the ground states. The coefficient 18 
means that the number of ground states of the Ising model is 18. Due to the two to one 
correspondence noted in the last section, the number c of the toric generators is 9. 
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So far we have calculated the number c of the toric generators by computing the partition 
function of the Ising model. In the language of the tiling problem, c is the number of tilings 
of the torus by three types of rhombi. In fact, however, this is not the only information we 
can extract from the calcuration of the partition function; we are able to count the numbers 
of X- y- and z- type rhombuses contained in each tiling. For this purpose, we consider which 
rhombus corresponds to each plaquet in Figure 11. For example, the top left plaquet (1-1 
element) in Figure 11 does not corresponds to any combination of rhombuses since the three 
spins of each triangle have the same spins, while the plaquet right to this one (1-2 element) 
corresponds to a combination of the upper half of the x-tjpe rhombus and the right half 
of the y-tjpe rhombus. In this way, we obtain combinations of rhombuses as depicted in 
Figure 12. 



nothing 



nothing 



nothing 



nothing 



nothing 



nothing 



Figure 12: Combinations of rhombuses corresponding to each plaquet in Figure 11. 



Now we assign weights that count the number of rhombuses to each element in Figure 
12. For example, we assign zero to the top left element and y^^ry to the right to this one. 
These data are summarized in in a matrix as: 



To 








^/xy 





Vxx 


\ 









z 


v^ 






\/xy 


z 


Vyy 







v 







V^ 





/ 



(4.26) 



Now we introduce a function Z^xn obtained by replacing T2 in Znxn with T2. Then, for 
example, we obtain 

1 - 



^2x2 



(4.27) 
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where we have omitted terms which do not correspond to the ground states. The first term 
x'^, for example, imphes that there is a tihng with four x-tjpc rhombuses, and the last 
term 2x'^y'^ implies that there are two tilings with two x-tjpe rhombuses and two y-tjpe 
rhombuses. Note that the numbers of x- y- and z- type rhombuses correspond to the three 
entries of the toric generators hx in the equation (2.30), and hence we are able to draw 
toric diagrams. In Figure 13, the numbers of toric generators on each vertex are depicted 
for n — 3,4, 5. 



3t 

si 



\ 3 
• •. 

21 



1 3 



3 1 



4 ♦ 
6 ♦ 
4 • 



\ 6 

124 

• • 

124 124 



6 



5 * 

10 
10 
5 

1 



605 



10 



10 



• • • 

1905 1905 



605 1905 605 

■■• • • •■ 

5 10 10 5 



Figure 13: The numbers of toric generators on each vertex. 



So far we have considered only the models of the form C^/Zn x Z„, but it is straight- 
forward to generalize the analysis to the models of the form C^/ x Z^ The partition 
function for the model ^nxm is given by 

7 _ ^(PP) I 7(pa) I ^(ap) I r7(aa) r,Q\ 

where 

ZSI = Tr(TiP))-, (4.29) 

Z^::l = Tr[(T(p))-i?„], (4.30) 

= Tr(Ti^))-, (4.31) 

Zi^l = Tr[{TjtYRn]. (4.32) 

The calculation of the numbers of rhombuses of each type is also straightforward. Table 1 
shows the number c of homogeneous coordinates for the model C^/ Z„ x for several n 
and m. 

The analyses is also applicable to the orbifold C'^/ Zn- The McKay quiver of the orbifold 
C^/Zn is obtained from that of the orbifold C^/Zn x Z„ by the identification [10], 

Ri,j = Ri-a,j+l- (4.33) 

Thus if we consider configurations of spins which satisfies this identification, we are able to 
compute the toric generators and the degeneracy. 
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1 


2 


3 


4 


5 


6 


1 


3 


5 


9 


17 


33 


65 


2 


5 


9 


17 


33 


65 


129 


3 


9 


17 


42 


113 


309 


860 


4 


17 


33 


113 


417 


1537 


5793 


5 


33 


65 


309 


1537 


7623 


38405 


6 


65 


129 


860 


5793 


38405 


263640 


7 


129 


257 


2445 


22081 


193849 


1849073 


8 


257 


513 


7073 


84609 


984577 


13034049 


9 


513 


1025 


20706 


326273 


5022543 


92168924 


10 


1025 


2049 


61097 


1265793 


25659185 


657904329 


11 


2049 


4097 


181245 


4933121 


131347909 


4735839173 



Table 1: The number c of homogeneous coordinates for the model C^/Z„ x Z. 
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